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Schramm-Loewner Evolution and isoheight lines of correlated landscapes 
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Real landscapes are usually characterized by long-range height-height correlations, which are 
quantified by the Hurst exponent H. We analyze the statistical properties of the isoheight lines for 
correlated landscapes oi H € [—1,1]. We show numerically that, for H < 0 the statistics of these 
lines is compatible with SLE and that established analytic results are recovered for H — —1 and 
H = Q. This result suggests that for negative H, in spite of the long-range nature of correlations, the 
statistics of isolines is fully encoded in a Brownian motion with a single parameter in the continuum 
limit. By contrast, for positive H we find that the one-dimensional time series encoding the isoheight 
lines is not Markovian and therefore not consistent with SLE. 
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We study isoheight lines of long-range correlated land¬ 
scapes. They are the paths of constant height in to¬ 
pography m, the equipotential lines on energy land¬ 
scapes [ 213 ], and the constant vorticity lines in turbu¬ 
lent vorticity fields [4]. Empirical and numerical studies 
of isoheight lines show that they are usually scale invari¬ 
ant and that their fractal dimension df depends 

on the long-range correlations, quantified by the Hurst 
exponent H |7]. 

One solid framework to study fractal curves in two- 
dimensions is the Schramm-Loewner Evolution (SLE) 
theory [S]. Accordingly, an SLE curve can be mapped 
onto a one-dimensional Brownian motion of a single pa¬ 
rameter in the continuum limit. Establishing that iso¬ 
height lines are SLE would allow us to generate an en¬ 
semble of such curves by simply solving a stochastic dif¬ 
ferential equation, without generating the entire land¬ 
scape, and to extend established results from SLE to 
isoheight lines. 

We analyze the zero isoheight lines of random surfaces 
of different H generated with free boundary conditions 
on a triangular lattice of size L^ x Ly with L^ > Ly. We 
consider chordal curves that are non-self-touching curves 
growing to infinity in the upper half-plane H, starting at 
the origin. From the Riemann mapping theorem, there 
is a unique conformal map gt that iteratively maps the 
complement of such a curve in the upper half-plane El 
back onto El. This map satisfies the Loewner differential 
equation, 

dgtjz) 2 

dt 

with 50 (x) = z and a real continuous function, called 

the driving function. If a curve is SLE^, then 

where Bt is a standard one-dimensional Brownian motion 


and K is the diffusion constant. We show that, while for 
H < 0 the statistics of is consistent with a Brownian 
motion, for H > 0, is not Markovian and therefore 
SLE cannot be established. 

We generate random landscapes on a triangular lattice, 
by assigning to each lattice site x = (x, y) a random 
height /i(x), imposing long-range correlations with the 
following spectrum, 

^(q)-|qr''= , (2) 

where = ‘^{H + 1) and H is the Hurst exponent. For 
that, we use the Fourier Filtering Method isHm where, 

h{x)=E (^yS'(q)u(q)) , (3) 

where J- denotes the inverse Fourier transform and 
the M(q) are independent complex Gaussian random 
variables of mean zero and unitary variance satisfying 
m(— q) = rt(q). With this scheme, one recovers uncorre¬ 
lated landscapes for H = — 1 and the discrete Gaussian 
Free Field (GFF) for iJ = 0 [H]. 

To obtain the zero isoheight lines, we extract the set 
of bonds on the dual lattice separating the sites of neg¬ 
ative height from the positive ones. We then focus on 
their accessible perimeter obtained by moving along the 
isoheight line and shortcutting distances equal to the lat¬ 
tice unit of the dual lattice, i.e., -x/S/S lattice units of the 
triangular lattice [2. Note that, using the formalism of 
ranked surfaces, one can show that, for H < 0, these lines 
correspond to the accessible perimeter of a correlated per¬ 
colation cluster at the percolation threshold [THIS]. The 
accessible perimeters of the zero isoheight lines on the 
triangular lattice in the case oi H = — 1 and H = 0 are 
analytically tractable and they have been proven to be 
SLE^ and SLE 4 , respectively [a Hang. 
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SLE and fractal dimension. As proven by Bef- 
fara [16], SLE^^ curves are fractals of a fractal dimension 
df that is related to the diffusion coefficient k by, 

df = min { 2 ,1 + 0 ■ (4) 

The fractal dimension d/ of the accessible perimeter of 
the isoheight lines for different values of H was numeri¬ 
cally estimated and even a conjecture was proposed for 
its dependence on H |3I7|. Using Eq. Q, this gives a 
first estimate for the expected values of k if the curves 
are SLE, see Table|lj To verify if SLE can be established 
we will compare the k values calculated from df with es¬ 
timates obtained with two indirect methods, the winding 
angle and the left-passage probability, and with the one 
obtained from the direct SLE mapping. 

Winding angle. The winding angle of SLEf^ curves 
follows a Gaussian distribution and the variance scales 
with K as, 

{e^)-{ef = al = b + -^HL,), (5) 

where & is a constant and Ly is the vertical lattice size 
[HiiniiiHj- To verify this relation for each path, we con¬ 
sider the discrete set of points Zi of the path. The winding 
angle 9i at each point Zi can be computed iteratively as 
9i+i — 9i + ai, where ai is the turning angle between two 
consecutive points Zi and on the path. Duplantier 
and Saleur computed the probability distribution of the 
winding angle for random curves using conformal invari¬ 
ance and Coulomb gas techniques m- k/A corresponds 
to the slope of ag against ln{Ly), see Fig. The esti¬ 
mates of K are displayed in Table |lj 

For values near H = 0 and H = 1, one has less preci¬ 
sion on the results, as the system is strongly influenced 
by finite-size effects, see e.g. Ref. |7|. The results we 
obtain from the winding angle measurement are, within 
error bars, in agreement with the previous estimates from 
the fractal dimension of the accessible perimeter of zero 
isoheight lines. Indeed, Eq. © gives insights into the 
conformal invariance of the problem [T]. Our results in 
Fig.[T] give some numerical indication that the accessible 
zero isoheight lines display conformal invariance, which 
is a prerequisite for SLE. 

Left-Passage Probability. As we simulate the 
curves in a bounded rectangular domain, we map confor¬ 
mally the isoheight lines into the upper half-plane, using 
an inverse Schwarz-Christoffel transformation HU, to ob¬ 
tain the chordal curve, which splits the system into two 
sides. For chordal SLE curves, Schramm has computed 
the probability Pk,{4 ') that a given point z = in the 
upper half-plane H is on the right-hand side of the curve 
m- This probability only depends on (j) and is given by 
Schramm’s formula 

Q- ■ 
( 6 ) 



FIG. 1. (color online) Rescaled variance of 

the winding angle for different Hurst exponents = 
-0.7,-0.55,-0.25,0,0.25,0.55,0.7. The black solid line rep- 

a-^—b 

resents the relation ^. = lnL„. In the inset, the rescaled 
probability distributions are plotted and compared to a nor¬ 
mal distribution for H = —0.7, 0,0.7. 

where F is the Gamma function and 2 F'i is the Gauss 
hypergeometric function [21j . This probability is known 
as left-passage probability. 

We define a set of sample points S' in El for which we 
measure the left-passage probability in order to compare 
it to the values predicted by Schramm’s formula (|^. To 
estimate k, we minimize the mean square deviation Q{k) 
between the computed and predicted probabilities, 

= ( 7 ) 

where (j) = arg(z), |S| is the cardinality of the set S, and 
P{z) the measured left-passage probability at z. The 
estimated value of k corresponds to the point where the 
minimum of Q{k) is observed. 

As shown in Fig. for FI < 0, the minimum of Q is 
less pronounced for higher values of H, as it is expected 
for functions of the form © with values of k increasing 
towards four. As summarized in Table |T| the estimated 
values of n obtained for negative LI are consistent with 
the ones predicted from d / and confirmed by the winding 
angle analysis. However, for H > 0, this is not the case. 
The obtained k values do not significantly depend on H 
and are consistently higher than the ones obtained from 
df. This result suggest that for iF > 0 the curves are not 
SLE. 

Direct SLE. To further test if the curves are SLE, 
one has to check that the statistics of the driving function 

is consistent with a one-dimensional Brownian motion 
with variance at. This can be done by solving Eq. 0 nu¬ 
merically. In order to do so, we use the so-called vertical 
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FIG. 2. (color online) Measured rescaled mean square devi¬ 
ation Q{k)/Q{K min) as a function of with Kmin the 

value of K where the minimum of (5 (k) is attained, for differ¬ 
ent Hurst exponents H = —0.7, —0.55, —0.25, 0, 0.25,0.55, 0.7. 
We chose 50^ points in the range [—0.025Li,,O.O 25 L 2 ,] x 
[0.15Z/y,0.25Ly] with Ly = 1024 and = SLy, which are 
then mapped to the upper half-plane through an inverse 
Schwarz-Christoffel transformation m- 


xlO'^ 



FIG. 3. (color online) Rescaled variance of 

the driving functions for different values of H = 
—0.7,-0.55,-0.25,0,0.25,0.55,0.7. In the inset, we present 
the rescaled probability distributions of the driving func¬ 
tions and compare them to a Gaussian distribution for H = 
-0.7,0,0.7. 


slit map algorithm, where one considers the driving func¬ 
tion to be constant over small time intervals 6t. Making 
this approximation, one can solve Eq. ([^ to obtain the 
following slit map equation [22j [23] , 


(c(t,r))( of the increments, with 

(S^t+rS^t) - {S^t+r){S^i) 


z{t,T) = 




,V) - (<56+x)2) m?) - 


(9) 


gt{z) =^t + \J[z- + 4:5t. ( 8 ) 

At t = 0, one considers the initial curve consisting of the 
points {zq = 0 ,..., z% = zn}, and sets the driving func¬ 
tion to be ^0 = 0. Then at each iteration f = 1,..., A^, we 
apply the conformal map gt^ to the remaining points Zj~^ 
of the curve, for j = j,..., iV to obtain the new mapped 
curve. One gets a new set of points z^y+i = for 

j = i,..., N— 1 shorter by one point, by mapping zl~^ to 
the real axis. For that, in Eq. ® we set 

and Sti = ti — ti-i = (lm{zj“^}) /4, where Re{} and 
lm{} are the real and imaginary parts, respectively. 

We extracted the driving function of all paths and 
computed the diffusion coefficient k, see values in Ta¬ 
ble |lj from the variance of the driving function and tested 
its Gaussian distribution at a fixed Loewner time t, see 
Fig-i We find a linear scaling of the variance of the 
driving function with the Loewner time. However, only 
for H < 0, the estimated values of k are consistent with 
the ones predicted from df. For H > 0, we obtain values 
of K significantly higher than the ones expected from df. 
In fact, the value of k even increases with H, instead of 
decreasing as expected. 

We also test the Markovian property of the driv¬ 
ing function by computing the auto-correlation c(r) = 


as shown in Fig. where S^t is the change in the driving 
function at time t. For H < 0, the correlation func¬ 
tion drops to zero after few time steps, as expected for 
a Brownian motion. For H > 0 we observe that the 
correlation function decays as a power law and thus the 
driving function is not Markovian, as it shows persistence 
in the increments. This explains why, though the vari¬ 
ance of the driving function is a linear function of time, 
the statistics of the isoheight lines for iL > 0 are not 
consistent with SLE. 

Conclusion. We numerically showed that the statis¬ 
tics of the accessible perimeter of the zero isoheight lines 
of long-range correlated landscapes are consistent with 
SLE only for —I < < 0. For this range, results for 

the fractal dimension, winding angle and direct SLE are 
in agreement within error bars, see Fig. This means 
that one can describe these curves with a Brownian mo¬ 
tion parameterized by a diffusivity k. In the two analytic 
limits H = —1 and H = 0, the accessible perimeters are 
SLE^/s and SLE 4 for iJ = — 1 and H — 0, respectively, 
results that we recover in our numerical analysis. To our 
knowledge, this is the first time that for an entire range 
of values of the Hurst exponent 7L, a family of curves 
coupled to random landscapes is shown to be consistent 
with SLE. This gives new insight in the field of fractional 
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FIG. 4. Auto-correlation function c(r) of the increments for 
different values oi H = -0.7,-0.55,-0.25,0,0.25,0.55,0.7, 
averaged over 50 time steps. The solid lines are power laws 
of exponent —0.4. 


TABLE I. Diffusion coefficient k computed from the fractal 
dimension Kfrac using data from [7] for H < 0 and obtained 
numerically for H > 0 using the yardstick method, from the 
winding angle Ke, from the left-passage probability klpp and 
the direct SLE method KdSLE, for the different values of the 
Hurst exponent H. 


H 

^frac 


no 


«:lpp 

«:dSLE 

-1 

2.76 

± 

0.16 

2.66 

± 

0.01 

2.69 

± 

0.08 

2.66 

± 

0.12 

-0.85 

2.80 

± 

0.24 

2.76 

± 

0.02 

2.80 

± 

0.07 

2.79 

± 

0.10 

-0.7 

2.96 

± 

0.24 

2.94 

± 

0.03 

2.95 

± 

0.11 

2.97 

± 

0.17 

-0.55 

3.20 

± 

0.24 

3.14 

± 

0.03 

3.13 

± 

0.15 

3.29 

± 

0.32 

-0.4 

3.35 

± 

0.31 

3.32 

± 

0.03 

3.27 

± 

0.17 

3.46 

± 

0.26 

-0.25 

3.64 

± 

0.28 

3.45 

± 

0.09 

3.40 

± 

0.21 

3.67 

± 

0.26 

-0.1 

3.90 

± 

0.18 

3.49 

± 

0.04 

3.50 

± 

0.24 

3.84 

± 

0.20 

0 

3.88 

± 

0.20 

3.44 

± 

0.18 

3.52 

± 

0.25 

3.98 

± 

0.19 

0.1 

3.44 

± 

0.40 

3.44 

± 

0.27 

3.56 

± 

0.26 

4.07 

± 

0.13 

0.25 

3.04 

± 

0.24 

3.07 

± 

0.07 

3.59 

± 

0.29 

4.26 

± 

0.13 

0.55 

1.92 

± 

0.24 

2.22 

± 

0.08 

3.62 

± 

0.36 

4.84 

± 

0.29 

0.7 

1.36 

± 

0.16 

1.65 

± 

0.16 

3.62 

± 

0.41 

5.35 

± 

0.25 

0.95 

0.56 

± 

0.16 

0.62 

± 

0.31 

3.63 

± 

0.50 

5.97 

± 

0.51 
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FIG. 5. (color online) Estimated diffusion coefficients k 
from the fractal dimension, the winding angle, the left-passage 
probability (Ipp), and the direct SLE methods (dSLE). The 
red crosses correspond to the rigorous results. 


Gaussian Fields in two dimensions [12], what might be 
helpful for understanding correlated landscapes from a 
theoretical point of view. One might also wonder if these 
isoheight lines can be related to some random walk pro¬ 
cess, as in the case of uncorrelated random landscapes 
and the Gaussian Free Field. For example, the isoheight 
lines for iL = — I and H — 0 are two specific cases of the 
overruled harmonic walker [24j . 

For positive H, we show that the driving function also 
scales linearly in time but it is not Markovian, a necessary 
condition to establish SLE. The numerical data suggest 


that the auto-correlation function scales as a power law. 
Gredidio et al. showed that if the driving function is a 
stochastic process with anomalous diffusion, the gener¬ 
ated tracers are anisotropic [25] . In the same spirit, it 
would be of interest to systematically study the statis¬ 
tics of random curves generated from a driving function 
with persistence and compare them to isoheight lines. 

This Letter also opens the possibility of applying SLE 
to the study of landscapes with negative Hurst exponent 
H. As we have shown, many systems can be considered 
from the point of view of a landscape, where isoheights 
are of relevance. Especially, it has been shown that zero 
vorticity isolines in two-dimensional turbulence are SLE 
|3|. There have been also attempts to extend this result 
to isolines in a generalized Navier-Stokes equation [25] to 
study the conformal invariance of a larger class of turbu¬ 
lence problems. It would be interesting to see if a relation 
between this problem and our results can be drawn for 
the accessible perimeters of these contour lines. 
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